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Abstract: Recent numerical studies of the 4D pure compact U(l) lattice gauge 
theory, I have participated in, are reviewed. We look for a possibility to construct an 
interesting nonperturbatively renormalizable continuum theory at the phase transi- 
tion between the confinement and Coulomb phases. First I describe the numerical 
evidence, obtained from calculation of bulk observables on spherical lattices, that 
the theory has a non- Gaussian fixed point. Further the gauge-ball spectrum in the 
confinement phase is presented and its universality confirmed. The unexpected re- 
sult is that, in addition to massive states, the theory contains a very light, possibly 
massless scalar gauge ball. I also summarize results of studies of the compact U(l) 
lattice theory with fermion and scalar matter fields and point out that at strong 
coupling it represents a model of dynamical fermion mass generation. 
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1 Introduction 



Quantum field theories used in the standard model and its supersymmetric exten- 
sions are either asymptotically free or so-called trivial theories. Both are defined in 
the vicinity of Gaussian fixed points and can be studied perturbatively. In various 
collaborations we address the question whether in four dimensions (4D) there exist 
quantum field theories which are not accessible to or anticipated by perturbation 
theory. I describe the results of our systematic numerical study of critical behavior 
in several compact U(l) gauge models on the lattice, both pure and with matter 
fields. Their common feature is confinement when the bare gauge coupling is strong. 
It has never been clear whether this is only a lattice artefact or whether a confining 
continuum U(l) gauge theory with confinement can be constructed. Though not 
yet conclusive, our results are promising and we hope to stimulate more theoretical 
attention. 

2 Renormalization of lattice field theories 

In numerical simulations of lattice field theories the concept of renormalizability 
requires the existence of critical behaviour somewhere in the bare coupling parameter 
space. Dimensionful observables, e.g. masses m s , mi, rri2, ■ ■ ■ , n-point functions, 
etc, can be calculated in the lattice constant units a as functions of couplings in 
the vicinity of a critical point or manifold. The dimensionless correlation lengths 
& = l/am,i — > oo when a critical point is approached. 

Renormalizability further requires the existence of "lines of constant physics" 
in the bare coupling parameter space, along which the observables scale, i.e. their 
dimensionless ratios r\ = mi/m s ,r2 = m2/m s , . . . , stay (approximately, in practice) 
constant. If such a line hits the critical point, one can construct the continuum limit 
and obtain the values of dimensionful observables by fixing one mass scale, m s , in 
physical units. Then a = l/£ s m s — > and = Tjm s . These results are usually 
universal, i.e. independent of the detailed choice of the coupling space and lattice 
structure, being governed by a limited number of fixed points. 

The lines of constant physics can approach a critical manifold but, before hitting 
it, leave the parameter space. Or they can enter phase transitions of weak first 
order. In both cases the correlation lengths grow, but stay finite. Some inherent 
cut-off A = £™ ax m s remains. If it is large with respect to the physical scale m s , i.e. 
if t^ ax » 1, the theory is renormalizable in a restricted sense. In this way e.g. the 
familiar trivial theories arise. 

3 Pure gauge theory on spherical lattices 

Recently we have reconsidered the oldest candidate for a non-Gaussian fixed point 
in the 4D lattice field theory, the phase transition between the confinement and the 
Coulomb phases in the pure compact U(l) gauge theory. We have used the extended 
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Wilson action in order to enlarge the possibility of finding a second order phase 
transition, 

S = -Y,w P [/3cos(Qp) +7cos(26 P )] . (f) 
p 

Here wp = 1 and 0p G [0, 2tt) is the plaquette angle, i.e. the argument of the 
product of U(l) link variables along a plaquette P. Taking Q P = a 2 gF^ u , where a 
is the lattice spacing, and (3 + 47 = 1/g 2 , one obtains for weak coupling g the usual 
continuum action S — ^ / d A xF 2 v . 

Earlier investigations performed as usual on toroidal lattices suggested that the 
phase transition, which is clearly of first order at positive 7, might be of second order 
at negative 7. Though a weak two-state signal is present there [|l|, it might be a 
finite size effect. It has been found that it disappears on lattices with sphere-like 
topology 0, [| [|, ||]. As this type of lattice does not change the universality class 
||. |7|], one can use it as well as the usual toroidal one. 

It has turned out that on spherical lattices the transition has properties typical 
for a second-order transition. This concerns, in particular, the dependence on the size 
of the lattice. Use of modern finite size scaling (FSS) analysis techniques, and large 
computer resources, allowed to determine on finite lattices the critical properties 
of the phase transition. A more detailed account of our work, as well as relevant 
references, can be found in Refs. [Q, |4]] . Here I list only the most important findings. 

The measurements have been performed at 7 = 0, —0.2, —0.5 on lattices of the 
volumina up to about 20 4 . The FSS behavior of the Fisher zero, specific heat, some 
cumulants and pseudocritical temperatures gave consistent results. The value of the 
correlation length critical exponent v has been found in the range v = 0.35 — 0.40. 

The most reliable measurement of v has been provided by the FSS analysis of the 
Fisher zero, i.e. of the first zero zq of the partition function in the complex plane of 
the coupling (3. The expected behavior with increasing volume V is Imzo oc V~ x l Dv . 
The joint fit to the data at all three 7 values gives v = 0.365(8). The scaling 
behavior of various pseudocritical temperatures, which have been determined from 
several other observables, is consistent with this value, as seen in Fig. |l|. 

The consistency of the data with the FSS theory suggests that the phase transition 
is of second order, as it implies growing correlation lengths. Unfortunately, it does 
not exclude that this growth stops on lattices substantially larger than those we 
have used. Very weak first order phase transition is still possible [|[ ||, which would 
mean that the cut-off cannot be made really infinite. The scaling behaviour could 
be governed by a fixed point "behind" the phase transition. 

In a very recent paper (which appeared after the Symposium) |I| Campos et al. 
suggest that this is what happens. Of course, also their data require an extrapolation 
to the infinite volume. For example latent heat, associated with the two-state signal 
present on toroidal lattices, decreases with volume. It can be extrapolated both 
to nonvanishing values and zero, depending on the ansatz. At some second-order 
transitions, two-state signals due to finite size effects are known to vanish extremely 
slowly with volume, mocking up a first order. It is most probably not possible to 
clarify the issue beyond any doubt by numerical methods. But the theory might be 



3 



1.17 




1 1 


i 1 i i i i 1 








*:?: 


jf t = -0.2 




1.16 
















I, ''.. 




1.15 




o : 


0(Re(z Q )) \ - "-•< 








& : 


3(c v ) "■,! "-• 








x : 










v : 


3(U 4 ) s 




1.14 












.005 .01 

i /I 1/0.365 



Figure 1: Consistency of the FSS of various pseudocritical temperatures with the 
value v = 0.365 of the correlation length exponent v at 7 = —0.2. The figure is from 
Ref. § 



interesting even if the cut-off could not be made infinite but only much larger than 
the physical scale, as in physical applications we are dealing with effective theories 
anyhow. 



4 Gauge-ball spectrum in the confinement phase 

In any case the scaling behaviour of as many observables as possible should be 
determined. Therefore we have investigated at 7 = —0.2 and —0.5 the spectrum of 



the theory in the confinement phase |TT], [L2| . It consists of massive "gauge balls", 
states analogous to glue balls in pure QCD. In particular, no massless photon is 
present in this phase. Also the string tension a has been estimated. 

The gauge-ball masses rrtj in various channels j of the cubic group have been 
measured at both 7 for different (3. Then their scaling behaviour in the form 

m^cM-Pr (2) 

(and similar for y/o) has been determined in each channel j individually. We have 
found two groups of masses with strikingly different scaling behaviour. Within the 
whole (3 range a large group of the gauge-ball masses scale with roughly the same 
exponents Uj close to the non-Gaussian value 0.365(8) found in 0, |], [|]. Though its 
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Figure 2: Values of Uj obtained at 7 = —0.2 in each gauge-ball channel separately. 
The double vertical line separates two groups with distinctly different Vj. The results 
to the right of the dashed vertical lines correspond to heavier states which are very 
difficult to measure. The figure is from Ref. ||T3| 



accuracy is not yet satisfactory, the exponent of \fa seems to be consistent with this 
value, too. However, in several channels getting contribution from the ++ gauge ball 
the values of Uj are approximately Gaussian, i.e. 1/2. This is shown in Fig. |2|. The 
values Pi are quite consistent with each other in all channels for each 7. Therefore, 
in the further analysis we have assumed the same value of P c in all channels. 

Assuming the same exponent Uj for each group, denoted u ng and v g for the non- 
Gaussian and Gaussian group, respectively, a joint fit was performed with these two 
exponents, common P c , and the individual amplitudes Cj as free parameters: 

m o = c o rV ^ f = n g>g (3) 
The resulting values of the exponents at 7 = —0.2 are 

u Qg = 0.367(14) 

v g = 0.51(3). (4) 

The preliminary results at 7 = —0.5 are consistent with these values [Q. 

From these results we conclude that the system has two mass scales which we 
denote by m ng and m g . In Fig. |3] we show their behaviour. Except ++ , the gauge- 
ball masses and presumably \fa scale proportional to m Qg . The ratios r,- = mj/m ng 
are thus constant. This holds for both 7 values. Thus, within the limits of numerical 
determination, we have found two lines of constant physics. We have further found 
that 77 are independent of 7, which indicates universality. This suggests that the 
scaling behaviour of the pure U(l) gauge theory belongs at different 7 to the same 
universality class governed by one non-Gaussian fixed point. The parameter 7 is 
irrelevant in this class as long as it is kept negative. 
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Figure 3: Masses of the ++ (lower points) and l + ~ glue balls in the confinement 
phase versus f3 at 7 = —0.2. The curves are fits (|3|) with the v values @. The figure 
is from Ref. JT3[ . 



Choosing m ng for a physical scale, we can consider a continuum theory (or a 
theory with large cut-off) with the spectrum given by the values of rj. They are 
given in UTTjl - As m g /m ng approaches zero, this theory would contain massless (or 
very light) scalar, which could possibly decouple. 

What kind of theory would it be? The procedure we have applied to the spectrum 
and a can be extended to all observables and, in general, to n-point functions. If 
scaling in corresponding powers of m ng is found, one can, in principle, construct the 
theory without having any continuum Lagrangian for it. Or it could be some non- 
polynomial Lagrangian in the fields corresponding to the states in the spectrum. As 
a would possibly be finite and nonvanishing in the units, it would be a confining 
theory. An interesting proposal for an approximate continuum Lagrangian has been 
made in Ref. and reported at this workshop by N. Sasakura. 

The pure U(l) lattice gauge theory with the Villain (periodic Gaussian) action 
belongs to the same universality class ||. Rigorous dual relationships valid for 
that action imply that also the following 4D models are governed by the same non- 
Gaussian fixed point: the Coulomb gas of monopole loops fll5| , the noncompact U(l) 
Higgs model at large negative squared bare mass (frozen 4D superconductor) J[B], [17] , 
and an effective string theory equivalent to this Higgs model JT8L |i~9j. 



5 Compact U(l) gauge theory with matter fields 

These findings raise once again the question, whether in strongly interacting 4D gauge 
field theories further non-Gaussian fixed points exist, which might possibly be of 
interest for theories beyond the standard model. The pursuit of this question requires 
an introduction of matter fields. Therefore we have investigated some extensions of 
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the pure compact U(l) gauge theory which might have interesting fixed points. We 
introduce fermion and scalar matter fields of unit charge, either each separately or 
both simultaneously. Because of space limits I give only a very brief overview of the 
results. 

The compact fermionic QED has some analogous promising properties in the 
quenched approximation f20| . In compact scalar QED we confirm the Gaussian 



behavior at the endpoint of the Higgs phase transition line ^TJ. In a theory with 
both scalar and fermion matter fields {xU<fi model), the chiral symmetry is broken, 
and the mass of unconfined composite fermions F = <$x is generated dynamically 



[2"2"|. It might be an explanation or alternative to the Higgs- Yukawa sector of the SM. 



In 2D and 3D [[24J] it belongs very probably to the universality class of the 2D and 
3D Gross-Neveu model, respectively, and is thus nonperturbatively renormalizable. 
In the 4D xU<ft model we demonstrate the existence of a tricritical point where the 
scaling behaviour is distinctly difrerent from the four-fermion theory, enhancing the 



chances for renormalizability p2L 25, 21 1. Here, apart from the massive fermion 



and Goldstone boson, the spectrum contains also a massive ++ gauge ball. The 
particular role of this state both in pure U(l) and in the xU<p model is remarkable 
and needs a theoretical explanation. 

I thank J. Cox, W. Franzki, C. B. Lang, T. Neuhaus, A. Seyfried, and P.W. 
Stephenson for collaboration on the topics I have described, and J. Ambjorn, I. 
Campos, and H. Pfeiffer for discussions. The computations have been performed 
on the computers of HLRZ Julich, RWTH Aachen, and KFU Graz. Our work was 
supported by Deutsches BMBF and DFG. 
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